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Abstract 

We consider a model of quantum field theory with higher derivatives for a 
spinor field with (<fip) 2 selfinteraction. With the help of the Bethe-Salpeter 
equation we study the problem of the two particle bound states in the " chain" 
approximation. The existence of a scalar bound state is established. 
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1 Introduction 



One of the basic problems in Quantum Field Theory (QFT) is related to the 
presence of divergences in the S-niatrix elements. The reason for their ap- 
pearance is well known: the propagators do not tend sufficiently quick to zero 
at large momenta. In standard QFT where one deals with equations of mo- 
tion for Fermi (Boze) fields which are of first (second) order in the derivatives, 
the only solution of this problem is given by the renormalization method. 
The core of this method consists in the following - due to the arbitrari- 
ness in the S'-matrix fl|], one is able to add to the Lagrangean counterterms, 
which compensate the divergences. But it is well known that this method is 
not universal - there exist nonrenormalizable theories (which contain infinite 
number of different types divergent Feynman diagrams). An example is the 
theory of the weak interactions, proposed by Fermi. Its nonrenormalizabil- 
ity has led after long investigations to the creation of the Weinberg-Sallam 
theory. 

One possible way to overcome the problem with divergences in QFT is to 
reject the restriction about the order of derivatives in the Lagrangean (or in 
the equation of motion) i.e., to consider Higher Derivative QFT (HDQFT). 
The main advantage of HD theories is that propagators have better behaviour 
at large momenta. That means that a theory containing a certain number 
of divergences in the standard case, in the HD case become (partially or 
completely) free from them. 

Higher derivatives (HD) are not something new in field theory. It is 
known since a long time that HD could be used for regularization; different 
reguralarizations are equivalent to the usage of HD ||. Regularization 
with higher derivatives is convenient to use in the analysis of Yang-Mills fields 
||]. HD have been used in classical electrodynamics || for elimination of the 
problem with the infinite field-energy of the electron. 

In the present paper we deal with a HDQFT for a spinor field with scalar 
selfinteraction. The equation of motion for this field is of third order in the 
derivatives. The coupling constant is dimensionless (with respect to the term 
with highest derivative), which is peculiar to the Yang-Mills models. An im- 
portant role in the model plays an additional mass-parameter with large (but 
finite) value. This parameter guarantees that the action is dimensionless and 
spinor fields have the standard dimension i.e. 3/2. As we shall see below, the 
propagator in our model coincides with a spinor propagator regularized by 
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Pauli-Villars, and the additional parameter plays the role of the regulariza- 
tion mass. That is why we suppose that this parameter is much larger than 
the mass of the physical field presented in the theory. After quantization, 
this parameter represents the mass of the spinor fields with nonpositive norm 
- the so called HD "ghosts" 0. 

In the nonrelativistic limit the equation of motion coincides with the non- 
linear Schrodinger equation. The basic hypothesis in our work is that (at least 
for low energies) the equation of motion has solutions of the Schrodinger type, 
i.e., two-particle bound state. The main goal of the present article is to verify 
this hypothesis. In order to do that, we use an analogue of the Bethe-Salpeter 
equation |?J. 



2 The model 



In the present paper we discuss a model of QFT with third order (in the 
derivatives) Lagrangean L ~ (p(D — M 2 )(d fl , y fl — m)(p. However, for the 
derivation of the commutation relations and Green functions it is more con- 
venient to work with the following (equivalent to the initial) second order 
Lagrangean (with an additional spinor field ip): 



VM 2 - m 2 v ^ r " ' ry 2 
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where Ad^B = Ad^B—d^AB. This Lagrangean contains an additional spinor 
field ip(x), and M is the arbitrary but finite mass parameter mentioned above. 
The coupling constant g is dimensionless and we shall assume it has arbitrary 
value and sign. 

The equations of motion for both fields (p and ip are: 

\JM 2 — m? M z — m z 

(d 4- m 2 ) 

' > = (irfdp + mW (2) 



2 — m 2 



After exclusion of the auxiliary field ip from (|2]) we obtain a covariant equation 
of motion for the physical field only: 

(n + M 2 )(iYd^-m)cp + 2g(v<p)cp = Q (3) 
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Let us consider its nonrelativistic limit. For this, it is convenient to rewrite 
(§) in the form: 

/ □ _|_ m 2 \ 2q 

' +1 {i-fd» - m)ip(x) + — 2 J <p<p)<p(x) = o 



\M 2 — m 2 J v ' " M 2 — m 

In momentum representation the first term has the form: 

— p2 _|_ ^2 \ 

M^r + iJ(p- *0v(p) 

Because M 2 >> m 2 , we see, that for low momenta the term ~f 2 +m , is ex- 
tremely small. Consequently, for low energies the behaviour of (|3|) is deter- 
mined by the equation: 

{irf% - m)<p + M 2 _ m 2 ^)v = (4) 

This is a non-linear Dirac equation. Following the standard procedure, it is 
easy to find its nonrelativistic limit: 

Here ip' is the large component of the spinor field (p with mass m. The equa- 
tion (|^) is equivalent to the two-particle Schrodinger equation with contact 
potential. The latter has the form: 

d 

i—^(x x , x 2 ; t) = H{x x - x 2 )^(x 1 ,x 2 ; t) (6) 



where: 



1 2 

H = -— E V? + X5(x 1 - x 2 ) =H + Him 



Equation @ has a unique solution when the coupling constant A has negative 
sign JJ. This suggests, that at least in the domain of low energies, the 
equation (Q) could describe bound states of two relativistic particles. 
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3 Quantization of the model 



The quantization of the model is made in the interaction representation. This 
procedure consists in the replacement of the field functions with operators. 
We assume all operators are normal ordered, and for simplicity we omit the 
sign for normal ordering. The free equation of motion for the field (p is: 

(□ + M 2 ) 

lt^W-mV = (7) 

The anticomutator T a p(x) = {<p a (x), <^g(0)}, obeys the following initial con- 
ditions (see Appendix A): 

F af 3(x)\ xo=0 = ; d T a p(x)\ Xo=0 = ; d^T a/3 (x)\ Xo= o = (M 2 - m 2 )^5 3 (x) 

(8) 

It is obvious that T(x) satisfy (^). The solution of (0) with initial conditions 
©is: 

r Q/3 (x) = -i{i^d^ + m) af3 (D m (x) - D M (x)) (9) 

where D m (x) is the Pauli- Jordan function for mass m. 

The propagator of the field ip, here denoted by S^^x) = (T((p a (x)ip~p(0))) , 
is the solution of the equation: 

(D + M 2 ) 

-{i-fdp - m)iS\x) = -8{x) (10) 



M 2 — m? 



In momentum representation the solution of ( |T0| ) is: (/c M 7 M +m) 
and the propagator is: 



l i 



m 2 ~k 2 M 2 -k 2 



S^{x) = -i(i-fd» + mU(D c m (x) - D c M (x)) (11) 

where: 

(27r) 4 J m z — k z — ie 

is the causal Green function of the scalar field. The advantage of the usage 
of HD is obvious from the form (|TTJ) of the S c : the propagators in our model 
coincide with the ones in the standard theory of a spinor field regularized by 
Pauli- Villars. Now, it is easy to convince ourselves that in the present theory 
there are no divergences connected to the Feynman diagrams. 
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The solution of (|7|) could be written as: 

+7^p / d*k{e^{N + u^{k)at{k) + N^{k)ct{k)) 
+e- ikx (N + u2-(k)a;(k) + N.v r a ~ (k)c; (k))} (12) 

where (p m (x) is the usual spinor field satisfying equation: (i r y tJ, d fl —m)ip m ( y x) = 
0. The mode-vectors u(k) and v(k) obey the equations: 

(=Ffc"7M - M)u r ^{k) = 
( T ^7 M + M)t» r ' ± (A;) = 

and the usual normalising conditions 

u r ^{k)u s ' T {k) = = v r ^{k)v s ^{k) 



L M a ' ( fc ) = ^ 



2ft 

E^W^W = (fc " 7At 2 t M)a " ( 13 ) 



where ft = \Jk 2 + M 2 and JV± = ^gp. The creat ion and annihilation 
operators a ± and c ± in ([12]) , obey the following anticommutation relations: 

{a~(k), a^p) } = {a~{k)i a+ (p) } = -<5 rs 5 (A; - p) 
{c r -(fc),c^(p)} = c+(p)} = -M(fc-p) (14) 

According to eq. (|i~4"|) these operators represent "ghosts" (with mass M) in 
the model. The physical subspace is given by vectors | ph), which satisfy the 
following additional conditions: 



ph) = a | ph) = 
ph) = c~| ph) = 0, 



5 



4 The bound state equation 



In 1951 Bethe and Salpeter derived an relativistic integral equation for the 
two-particle bound states J/J. With the help of this equation, some rela- 
tivistic corrections to the fine and the superfine structure of the light atoms 
have been calculated 0, [§]. The kernel of this integral equation consists 
of all two-particle irreducible Feynman diagrams, but in the case of a small 
coupling constant one can consider the kernel restricted to its first term 
(the so called "ladder" approximation). Applying their approach, we obtain 
an analogue equation for the interaction, considered in the paper. Here, this 
approximation leads to extreme simplification - due to the presence of a 
5-f unctions, the integral equation is reduced to a set of 16 linear algebraic 
equations. In this approximation, the kernel corresponds to contact interac- 
tion of two particles, and could be written as: 

K(xi, 2 2 , x 3 , x 4 ) = ^S(x 1 - x 3 )5(x 2 - x 4 )<5(2i - x 2 ) (15) 

Hence, the analogue of the Bethe-Salpeter equation is: 

XPifc (p; x 5, x s )dx 3 dx 4 dx 5 dx 6 (16) 
where xip] x ii x 2) is the relativistic two-particle wave function, defined as: 

X a p(p;xi,x 2 ) = (0 | T(v9 Q (a;i)v? / 3(a; 2 )) | p) (17) 
Here p is the 4-momenta of the system, and 

C aia2/ 3 1/ 3 2 (2i, X 2 , X 3 , X 4 ) = S^fc ( x l — x s)^aaP2 ( X2 ~ 

is the product of two free Green functions ( [TT]) corresponding to the two 
particles. Using the transformation properties of the wave function: 

X(p; x i, x 2) = e ipa x(p; x i + a, x 2 + a) 

(where a is an arbitrary 4-vector) and introducing absolute and relative co- 
ordinates X and x resp.: 

Xi + x 2 
2 

X = X\ — x 2 
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we obtain: x{P'i x ii x 2) = e x( x )- m the new co-ordinates, instead of 
eg . ( fltf ) we have: 



^ _ m2) / ^(^(P-^XftA^ = 0) 

(18) 

One can show, that in the nonrelativistic limit eq. ([Tj|) coincides with 
the Schrodinger equation with ^-function potential. Indeed, for low energies 
(see (0)) in the momentum representation one can rewrite flTE|) as: 



f| + q-m) (|-g-mj XPxpM = ~ m 2 ^ia 2 ( x = 0) 

\ / ai/3i \ / a 2 f3 2 

where g is the relative momentum. Acting on both sides of the above equation 
with the projectors A + (q)A~(q), where A ± (q) = \ + ^ we obtain, 

after simple transformations: 

| - E® + s - tf) (| - E(q) - e - tf) x + '~(^ g) = ~ M2 ^ m2 X + '> = 0) 
Here is the energy in the rest frame (p = 0), e is the relative energy, E(q) = 



(q) 2 + m 2 and x + ' (£> <?) = A + (g)A (g*)x( £ 5 ?*) (resp. for the function x(0))- 
The latter equation could be integrated with respect to e. The result is: 

(E - 2E(q)) J x + '-(e,q) de = ^—^ x +'-( x = 0) 

In the nonrelativistic limit the projectors could be removed from the above 
(scalar) equation and we obtain exactly the Schrodinger equation in the mo- 
mentum representation: 

Q 2 ~ 9 
(E - 2m)x(q) = —XW + jri o^ x = °) 



Let us return to the equation (|TE|). It is an algebraic system connecting 
the wave function, taken in different points (x and 0). Consequently, if we 
set x = we obtain the following set of algebraic equations: 

Xa ia2 {x = 0) = J r ai l3 1 a 2 /3 2 (p)X(3 1 f} 2 ( x = °) ( 19 ) 

where: 

^..Sim&IP) = (2 T )4(A/2 - , n 2) J ^.MffoftlP - k ) d ' k 

S%k) = -i^ + m )( 2 _\ 2 _. £ - M2 _\ 2 _ \ (20) 



The calculation of the integral (p0|) is long, but it is not connected to any 
principle difficulties. This integral is similar to the vacuum polarisation dia- 
gram in the Quantum Electrodynamics. Some details about the calculation 
of such integrals one can see in |J • 

The final result can be written in the following form: 

T = <V7^ ® l v + K i 1 ® 7 M + 7 M ® 1) + d 1 <g> I (21) 

where the functions a Mi/ , and d in the rest frame p = (E,0) are given in 
Appendix B. We are interested only in bound states between particles with 
mass m. This means, that the values of the energy must be in the interval 
< E 2 < Am 2 . It is easy to check, that in this interval a^, b^, d are analytic 
functions of the energy (the latter is not obvious only for the point E 2 = 0). 

If one writes the wave function x as column with 16 components, then 
the condition for compatibility of the set of algebraic equations ([19]) 

det(F l6xW - l 16xl6 ) = (22) 

determines whether there exist bound states or not. 

5 Decomposition into Lorentz invariants 

For the subsequent analysis it is instructive to decompose the two-particle 
wave function into Lorentz invariant quantities. First, we start with the 

observation, (see eq.(|TT|)) that the quantity x{x) = x( x ) (^ T ) > wnere C 
is the charge conjugation matrix, transforms as a tensor under bispinorial 
representation r(|, 0) © r(0, |) of the Lorentz group, i.e.: 

x'(x')=S- 1 x(x)S , x'^ = A^x„ (23) 

This representation is reducible, and it decomposes into a direct sum of 
irreducible representations: a scalar (s), a pseudoscalar (p), a 4- vector (V a ), 
an axial vector (A a ) and an asymmetric 2-rank tensor (S Q/ g). Using the 
connection S~ lr y^S = A%j u , x(x) could be written as: 

X = si + jrf 5 + V al a + A a7 V + (24) 

where a a/3 = \{ r ) OLr ) 13 — 7^7"). After simple transformations, from ( |2"T| ) and 
(|19"D, one can obtain the equation for x'- 

a^-fxi* + b» (xY - Yx) + (1 - d)x = (25) 
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where we use that C 1 ^f fl C = — 7J. Finally, after standard 7-matrix algebra, 
from (^5|) we obtain the following equations for the Lorentz invariants: 

(a£-d+l)s = 

(a£ + d-l)p = 2b^A p 

2a°V p - (a p + d-l)V a = 2b p iW a (26) 

2a°A p -{a p -d + l)A (T = 2Wp 

(a£-d+ l)tZ ap - 2{a aa iY, p a - a pa iY? a ) = 2{b a V p - b p V a ) 

Hence, we see that in terms of fields s,p, V, A, £ the matrix JF 16xl6 could be 
written as: 

^16x16 = Fixi(s) © F 5x5 (p, A) © ^i xio(V, S) 

Therefore, the condition for existing of bound states (p2|) decomposes to tree 
independent conditions: 

T{s) -1 = 0, det(.F(p, A) - 1) = , det{T{V, S) - I) = (27) 

6 Study of the scalar bound state 

Let us concentrate our attention on the conditions (p7|). It is clear that the 
condition for the scalar bound state is the simplest one, because the scalar 
field does not mix with the other fields. In this section we solve the condition 
for the scalar bound state, which (following (H)) is: 

a £-d+l = (28) 

This equation is too complicated to be solved in the general case. But since 
we have a small parameter - E 2 /M 2 - we can simplify eq.(^). For this 
purpose, we use the approximation of the functions a^ u and d given into 
Appendix B. Then eq.(|2"8"D is: 

1 + 9 M2 9 m2 \( E2 3)1 (—) 

16tt 2 M 2 - m 2 16tt 2 M 2 - m 2 \ \2m 2 J H \ m 2 ) 
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An interesting feature of this equation is that it contains two types of terms: 
the first two are of the order of 1, while the rest are ~ m 2 /M 2 and therefore 
are much smaller than the first two. This means that has solution if the 
large and small terms are equal to zero independently i.e. ( [29]) decomposes 
in two different conditions for the large and small parts. These two equations 
are: 

(30) 

M — - 1 




E 2 


-s) 


2m 2 




4) 


Rm 


m z 1 


~E 



E 
R~ 



(31) 



Remark : The separation of the eq. (p9|) in two equations is possible, if the 
large part is determined up to term of order m 2 /M 2 . This arbitrary term 
contributes in the equation for the small parts. In order to do that we 
represent t he coupling constant as: g = g' (l + Agj^J. From (j3~0D we have: 
g' = — 167T 2 . It is easy to verify that ( pT|) has a unique solution which depends 
on M 2 and Ag. Indeed, the LHS of ([H]) is a monotonously increasing function 
of E, and at the ends of the interval < E 2 < Am 2 it has the values: 



Ag + 3-3ln(Ml) E 2 = 
Ag-l-ln (K) E 2 = Am 2 



LHS§T{) = 

Consequently, the model considered here admits the existence of the scalar 
bound state with mass (in the rest frame) E 2 (Ag, M 2 ), which is the solution 
of eq. (|31"D . This bound state exists for each M 2 : M 2 » m 2 , when the 
coupling constant is negative and has the value: 

(77l 2 \ 
1 + A9 W 2 ) 

where: In (jf^) + 1 < Ag < 3ln (^^pj — 3. This result is also in correspon- 
dence with the solution of the non-linear Schrodinger equation, discussed 
above. 

In conclusion, we want to mark the possibility for solutions of (p7|), con- 
nected to the other fields. It is easy to verify from eq.(^), that in the rest 
frame, the matrices ^(p, A) and FiV, E), has only symmetrical placed nondi- 
agonal terms (besides the main diagonal) which are ~ m 2 /M 2 . Consequently, 
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their contribution to the determinant is ~ m 4 /M 4 or smaller. Because of this, 
we can consider these nondiagonal elements in these matrices as zeros. Fi- 
nally we obtain equations of the type for the four fields in consideration. 
To solve these equations we proceed in analogy with ( [30]) and ([H]). In order 
to satisfy the (corresponding) condition for the large parts, we have to choose 
the constant g' to be: 

g' p = 16tt 2 , g' v = 32vr 2 , g' A = -32vr 2 

The root, corresponding to £ field is not achievable, because the large com- 
ponent in this case is zero. It is important to note, that equations for small 
parts, corresponding to temporal and spatial components of the vector and 
axial- vector fields are different. 



Appendix A 



First, we introduce two auxiliary functions A a/3 (x) = {ip a (x), ^(O)} and 
S af 3(x) = {(p a (x),ij)p(0)}. 

From the requirement that fields, separated with space-like interval are 
causal independent, it follows, that: 

r<xp(%) U°=o = A a/3 (rr) 1^0=0 = S a p(x)\ x o =0 = (32) 

The procedure of canonical quantization gives: 

{n<p a (x),<pf}(0)} Uo =0 = {7^(20,^(0)} 1*0=0 = iS 3 (x)5 a p (33) 

where 7r is canonical conjugated momentum for corresponding field. Substi- 
tuting in (^) the explicit form of conjugated momenta, we have: 

d T a p (x) \ x o =0 = d A af3 (x) 1*0=0 = (34) 

and 

d S al3 (x) Uo=o = -zVM 2 - mH\x)5 aP (35) 
Let us consider eq.(|]) in the free case. One can show, that: 

Finally, from the latter equation and from eq. fl3"5] ) we obtain: 

d 2 T a p(x)\ x » =0 = (M 2 - m 2 ) 7 y 3 (o;) (36) 
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Appendix B 



In the rest frame p = (E, 0), functions a Mi/ , and d in (|2l|) are: 

(M 2 - m 2 ) 2 



9 I fiV 

16tt 2 E 2 (M 2 -m 2 ) \ 2 

(M 2 - m 2 ) 

+ (3E 4 — 3E 2 (M 2 + m 2 ) + (M 2 — m 2 ^ 



QE 2 



In 



M 2 



3E 2 



R 3 arth 



R 



£ n , ( E 

+ -R m arctg ^— 



M 2 + m 2 - E 2 

+ — R M arctg 



E 



R 



M 



+9^09u0 



2\2 



: (M 2 -m 1 ) 



2/^2 , 2^ o/^2 2^(^-™ 2 )- ^ 



+ (3E 2 (M 2 + m 2 ) - 2(AT - m 2 ) 2 ) 
1 



6-E 2 \ m 



3E 2 
E 



(E + E 2 (M 2 + m 2 ) - 2(M 2 - m 2 ) 2 )R arth 



R 



^{2m 2 + E 2 )R m arctg ( — 



M 2 + m 2 - E 2 

E / E 

— (2M 2 + E 2 )R M arctg 



M • 



d 



2m 2 g 



16n 2 E 2 (M 2 - m 2 ) 



(M 2 -m 2 ) ln (M 2 



R arth 



R 



M 2 + m 2 -E 2 , 



—ER m arctg 

E 
2m 



E 

Rm 



— ERMarctg 



E 



R 



where: 

R = JW 



M 2 



m 



2)2 _ 4M 2 m 2 R m = V4m 2 - E 2 R M = y/AM 2 - E 2 



If we recall that E 2 /M 2 is an extremely small parameter, then the above 
functions could be rewritten in the form: 
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"in/ 



16tt 2 (M 2 - m 2 ) 



9iiv 



M 2 m 2 



-In 



' M 2 



5 
6 
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) 



1_\ & 
18/ m? 



) 




lf arctg lid 

If E?\Rrn 

3 \ + ri 2 ~)~E 



2 2 

+9no9uom - 





d = 



16tt 2 (M 2 - m 2 ) 
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